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We present measurements of branching fractions and CP-violating asymmetries for neutral B 
meson decays to two-body final states of charged pions and kaons based on a sample of about 
88 million T(AS) — > BB decays. From a time-independent fit we measure the charge-averaged 
branching fractions B(B° -> tt+tt") = (4.7±0.6±0.2) x 10~ 6 , B(B° -c K + n~) = (17.9±0.9±0.7) x 
10~ 6 , and the direct CP- violating charge asymmetry An* = -0.102+0.050+0.016 [-0.188, -0.016], 
where the ranges in square brackets indicate the 90% confidence intervals. From a time-dependent 
fit we measure the B° -> tt+tv' CP-violating parameters = 0.02 ± 0.34 ± 0.05 [-0.54, +0.58] 
and CU = -0.30 ±0.25 ±0.04 [-0.72, +0.12]. 
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Recent measurements of the CP-violating asymmetry 
parameter sin2/3 reported by the BABAR Q and Belle Q 
Collaborations established CP violation in neutral B de- 
cays. These results are consistent with the Standard 
Model (SM) expectation based on indirect constraints 
on the magnitudes of the elements of the Cabibbo- 
Kobayashi-Maskawa Q quark-mixing matrix. However, 
a full test of the CP violation mechanism in the SM, 
through a single complex phase in the CKM matrix, will 
require additional direct constraints on the angles (a, /?, 
and 7) of the Unitarity Triangle 0. 

The time-dependent CP-violating asymmetry in the 
decay B° — > tt + it~ is related to the angle a, and ratios of 
branching fractions for various tttt and Ktt decay modes 
are sensitive to the angle 7. In this Letter we present 
results for branching fractions and CP-violating asym- 
metries in B a — > 7r + 7r~, K + ir~, and K + K~ decays Q 
using a sample of 87.9 ± 1.0 million BB pairs. A detailed 
description of the BABAR detector is presented in Ref. Q , 
and more details on the analysis technique are given in 
Refs. |7j, which describe our previous measurements of 
these quantities. Other measurements of the branching 
fractions and the charge asymmetry in B° — > K + ir~ have 
been performed by the CLEO and Belle Collaborations 
||. More recently, the Belle Collaboration reported a 
measurement of the time-dependent CP asymmetry in 

B° -» 7T+ 7T- B- 

We reconstruct a sample of neutral B mesons (P re c) 
decaying to the h + h!~ final state, where h and h' re- 
fer to 7r or K. Signal yields are determined with a 
maximum likelihood fit including kinematic, topological, 
and particle identification information. For the K^ir^ 
components, the yield is parameterized as N K ^ n ± = 
Nk-k (1 i Ak"7t) /2, where Nk-k is the total yield and 
A K « = (N K - n+ ~ N K+ „-)/{N K -„+ + N K+n -) is the 
CP-violating charge asymmetry. The asymmetry arises 
from interference between the b — ► s pen guin and b — > u 
tree amplitudes, and is predicted [lfj|ll( to be less than 
20% in the Standard Model. However, a larger asymme- 
try could be induced by new particles, such as charged 
Higgs bosons or supersymmetric particles, contributing 
to the penguin amplitude. 

In order to extract the CP asymmetry parameters in 
P° — > 7r + 7r~ decay, we examine each event in the P rec 
sample to determine whether the second B meson (Ptag) 
decayed as a B° or B° (flavor tag) and reconstruct the 
difference At between the proper decay times of the P rec 

?tag 

when h + h'~ = ir" 



parameters S-mr and CW are defined as 



and Ptag decays. The decay rate distribution /+ (/_) 
~7r~ and Ptag = P° (P°) is given by 



e -|At|/r 

/±(Ai) = — [1 ± S^sin(Am d At) 



4r 



T CWcos(Am d Ai)], (1) 



where r is the mean P° lifetime and Ato^ is the mix- 
ing frequency due to the eigenstate mass difference. The 



2ImX , _ 1 - |A| Z 
and Ctttt 



1 ■ 



i + iAr 



(2) 



and vanish in the absence of CP violation. If the decay 
proceeds purely through the b — > u tree amplitude, the 
complex parameter A is given by 



A(P 



v t iv u \ (V u * d V uh 



V V* \ V V* 



(3) 



In this case C^ = and — sin2a, where a = 
&rg[-V td V t l/V ud V* b }. In general, the b -> d penguin 
amplitude modifies both the magnitude and phase of A, 
so that Ctttt 7^ and S^n — yl — C% n sin 2a c ff, where 
a e s depends on the magnitudes and relative strong and 
weak phases of the tree and penguin amplitudes. Several 
approaches have been proposed to obtain information on 
a in the presence of penguins 0, . 

The event selection and P rec reconstruction used in 
this analysis are similar to those used in Ref. 0. 
Hadronic events are selected based on charged parti- 
cle multiplicity and event topology. Candidate P rcc de- 
cays are reconstructed from pairs of oppositely-charged 
tracks forming a good quality vertex, where the P rec four- 
momentum is calculated with the pion mass assumed for 
both tracks. 

Signal decays are identified kinematically using two 
variables, the difference AP between the center- 
of-mass (CM) energy of the P roc candidate and 
\fs/2, and the beam-energy substituted mass tties = 
y/{s/2 + pi ■ p B ) 2 /E'f - p|, where y/a is the total CM 
energy, and the P rec momentum pe and the four- 
momentum of the initial state (P^,pi) are defined in 
the laboratory frame. For signal decays AP and tties 
are Gaussian distributed with resolutions of 26 MeV and 
2.6MeV/c 2 , respectively. For ir + ir~ decays AP peaks 
near zero, while for decays with one or two kaons the 
AP peak position is parameterized as a function of the 
kaon momenta in the laboratory frame, with an average 
shift of —45 MeV and -91 MeV, respectively. The dis- 
tribution of toes peaks near the B mass. We require 
5.20 < toes < 5.29GeV/c 2 and |AP| < 0.15 GeV. 

Identification of h + h'~ tracks as pions or kaons is ac- 
complished with the Cherenkov angle measurement 9 C 
from a detector of internally reflected Cherenkov light. 
We construct charge-dependent double-Gaussian proba- 
bility density functions (PDFs) from the difference be- 
tween measured and expected values of 9 C for the pion or 
kaon hypothesis, normalized by the error ag a . The PDF 
parameters are measured in a sample of D* + — * P°7r + , 
D° — > K~ir + decays, reconstructed in data. The typical 
separation between pions and kaons varies from 8ae a at 
2GeV/c to 2.5cre c at 4GeV/c. 

We have studied potential backgrounds from other B 
decays and find them to be negligible. Backgrounds from 
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TABLE I: Average tagging efficiency e, average mistag frac- 
tion w, mistag fraction difference Aw = w(B°) — w(B°), and 
effective tagging efficiency Q for signal events in each tagging 
category. The quantities are measured in the -Bfl av sample. 



Category 


e 


(%) 


w 


(%) 


Aw 


(%) 


Q 




) 


Lepton 


9.1 


± 0.2 


3.3 


± 0.7 


-1.6 


± 1.3 


8.0 


± 


0.3 


Kaon I 


16.6 


± 0.2 


9.5 


± 0.7 


-2.8 


± 1.3 


10.7 


± 


0.4 


Kaon II 


19.8 


± 0.3 


20.6 


± 0.8 


-5.3 


± 1.3 


6.7 


± 


0.4 


Inclusive 


20.1 


± 0.3 


31.7 


± 0.9 


-2.6 


± 1.4 


2.7 


± 


0.3 


Untagged 


34.4 


± 0.5 
















Total Q 














28.4 


± 


0.7 



the process e + e~ — > qq (q — u,d,s,c) are suppressed by 
their topology. In the CM frame we define the angle 9s 
between the sphericity axis of the B candidate and the 
sphericity axis of the remaining particles in the event, and 
require |cos0s| < 0.8, which removes 83% of this back- 
ground. For these particles we also define a Fisher 
criminant J 7 ^ 0.53-0.60x^.p* + 1.27xX; i ]5 l * |cos(0*)| 2 
where p* is the momentum of particle i and 9* is the an- 
gle between its momentum and the B rcc thrust axis in the 
CM frame. The shapes of T for signal and background 
events are included as PDFs in the maximum likelihood 
fit. 

We use a multivariate technique to determine the 
flavor of the -B tag meson. Separate neural networks are 
trained to identify primary leptons, kaons, soft pions 
from D* decays, and high-momentum charged particles 
from B decays. Events are assigned to one of five mutu- 
ally exclusive tagging categories based on the estimated 
mistag probability and the source of the tagging infor- 
mation (Table Q). The quality of tagging is expressed in 
terms of the effective efficiency Q = ^2 k Ck(l — 2wfc) 2 , 
where ej, and u?k are the efficiencies and mistag prob- 
abilities, respectively, for events tagged in category k. 
Table [I] summarizes the tagging performance measured 
in a data sample -Bflav of fully reconstructed neutral B 
decays to D(*)~(tt + , p + , a^). The assumption of equal 
tagging efficiencies and mistag probabilities for signal 
7r + 7r~, K + ir~, and K + K~ decays is validated in a de- 
tailed Monte Carlo simulation. The background hypoth- 
esis have separate tagging efficiencies. 

The time difference At is obtained from the known 
boost of the e + e~ system and the measured distance be- 
tween the z positions of the B rcc and -Btag decay ver- 
tices. A detailed description of the algorithm is given 
in Ref. Q- We require \At\ < 20 ps and a At < 2.5 ps, 
where a At is the error on At. The resolution function for 
signal candidates is a sum of three Gaussians, identical 
to the one described in Ref. [lij . with parameters deter- 
mined from a fit to the Bflav sample (including events in 
all five tagging categories) . The background At distribu- 



tion is modeled as the sum of an exponential convolved 
with a Gaussian, with two additional Gaussians to ac- 
count for tails. Common parameters are used to describe 
the background shape for all tagging categories. We find 
that 96% of background events are described by an effec- 
tive lifetime of approximately 0.7 ps. 

We use an unbinned extended maximum likelihood fit 
to extract yields and CP parameters from the B rec sam- 
ple. The likelihood for candidate j tagged in category k 
is obtained by summing the product of event yield Ni, 
tagging efficiency e,^, and probability Vi t k over the eight 
possible signal and background hypotheses i (referring 
to 7r + 7r _ , K + n~, K~tt + , and K + K~ decays). The ex- 
tended likelihood function for category k is 



C k = exp 




The probabilities Vik are evaluated 



(4) 

as the product 



of PDFs for each of the independent variables Xj = 
{m ES , AE, T,e+,6~, At}, where 0+ and 9~ are the 
Cherenkov angles for the positively and negatively 
charged tracks. We use separate PDF parameters for 
0+ and 9~ to account for possible systematic differences. 
The total likelihood C is the product of likelihoods for 
each tagging category, and the free parameters are de- 
termined by maximizing the quantity In C. The fitted 
sample contains 26070 events. 

Signal yields are determined from a fit excluding tag- 
ging or At information. There are 16 free parameters, 
including signal and background yields (6 parameters); 
Ktt asymmetries (2); and parameters for the background 
shapes in m ES (1), AE (2), and T (5). Table [n] sum- 
marizes signal yields, total efficiencies, charge-averaged 
branching fractions, and Ak-k- In the efficiency calcula- 
tion we neglect possible effects due to final state radiation 
from the B rec decay products. The significance of Ak-k is 
2.0, where significance is defined as the square root of the 
change in — 21og£ when Ak-k is fixed to zero. These re- 
sults are consistent with our previous measurements 0, 
and with measurements from other experiments 0- For 
the decay B° — > K + K~ we measure a yield of only 1 ± 8 
events and so compute a Bayesian 90% confidence level 
(C.L.) upper limit on the branching fraction. Ref. Q 
gives a detailed description of the method used. 

The dominant sources of systematic error on the 
branching fraction measurements are from possible fit 
bias (determined in large samples of Monte Carlo simu- 
lated events), uncertainty in track and 9 C reconstruction 
efficiencies, and imperfect knowledge of the PDF shapes. 
The calculation of selection efficiencies using Monte Carlo 
simulated decays has been checked against control sam- 
ples in data and residual uncertainties are included in the 
systematic error on branching fractions. For Aktt the 
systematic error is dominated by the 9 C PDF shape and 
possible charge bias in track reconstruction. The total 
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FIG. 1: Distributions of rriES and AE for events enhanced in 
signal (a), (b) ty + tt~ and (c), (d) K^tt decays. Solid curves 
represent projections of the maximum likelihood fit, dashed 
curves represent qq and nn «-» Kn cross-feed background. 



systematic error is computed as the sum in quadrature 
of the individual uncertainties. 

Figure H shows distributions of mss and AE after se- 
lecting on probability ratios to enhance the signal purity. 
The solid curve in each plot represents the fit projection 
after correcting for the efficiency of the additional selec- 
tion (52% for 7T7T, 79% for Kir). 

The parameters and C^-n are determined from a 
second fit including tagging and At information, where 
the J3fl av sample is included to determine the signal pa- 
rameters describing tagging information and the At res- 
olution function. The At PDF for signal ir + ir~ decays is 
given by Eq. ^ modified to include Wk and Awk for each 
tagging category and convolved with the signal resolution 
function. We also take into account possible differences 
in reconstruction and tagging efficiencies between B° and 
B° mesons. The At PDF for signal K + ir~ events takes 
into account B°-B° mixing based on the charge of the 
kaon and the flavor of Stag- 

A total of 76 parameters are varied in the fit, includ- 
ing the values of S^n and C W7r (2); signal and background 
yields (5); Kir charge asymmetries (2); signal and back- 
ground tagging efficiencies (16) and efficiency asymme- 
tries (16); signal mistag fraction and mistag fraction dif- 
ferences (8); signal resolution function (9); and parame- 
ters for the background shapes in rriES (5), AE (2), T (5), 
and At (6). We assume zero events from B° — ► K + K~ 
decays and we fix Tgo and Arrid to their world average 
values ^5| . As a means of validating the analysis tech- 
nique, we determine r and Am ( i in the B rec sample and 




FIG. 2: Distributions of At for events enhanced in signal -kit 
decays with _B t a g tagged as (a) B° (N B o) or (b) B° (JVgo), and 
(c) the asymmetry [N B o — Ng ] / [N B o + Ngo] as a function 
of At. Solid curves represent projections of the maximum 
likelihood fit, dashed curves represent the sum of qq and Kn 
background events. 



find t = (1.56 ± 0.07) ps and Am d = (0.52 ± 0.05) ps" 1 . 
The combined fit to the B rcc and -Bfl av samples yields 

S„« = 0.02 ± 0.34 (stat) ± 0.05 (syst) [-0.54, +0.58] , 
C„ = -0.30 ± 0.25 (stat) ± 0.04 (syst) [-0.72, +0.12] , 

where the range in square brackets indicates the 90% 
C.L. interval taking into account the systematic errors. 
The correlation between S n7r and C nn is —10%. The 
signal yields determined in this fit are within 3% of the 
yields obtained from the time-independent fit. System- 
atic uncertainties on and CV^ are dominated by im- 
perfect knowledge of the PDF shapes and possible fit 
bias. We also evaluate multiplicative systematic errors 
(0.015), which are calculated at one standard deviation 
and summed in quadrature with the additive system- 
atic uncertainties. Figure [21 shows distributions of At 
for events with J5 tag tagged as B° or B°, and the asym- 
metry as a function of At for tagged events enhanced in 
signal 7T7r decays. 

In summary, we have presented updated measurements 
of branching fractions and CP-violating asymmetries in 
B Q — > 7T+7T+ A' + 7r+ and K + K~ decays. These results 
are consistent with, and supersede our previous measure- 
ments 0- We do not observe large mixing- induced or 
direct CP violation in the time-dependent asymmetry of 
B° —>■ ir + ir~ decays, as reported in 0- 
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TABLE II: Summary of results for total detection efficiencies, fitted signal yields Ns, charge-averaged branching fractions B, 
and Akw Branching fractions are calculated assuming equal rates for T(4£>) — > B B° and B + B~ . The upper limits for N K + K - 
and B(B° K+K~) correspond to the 90% C.L. 



Mode 


Efficiency (%) 


N s 


S(10" 6 ) 




A K * 90% C.L. 




38.0 ±0.8 


157 ± 19 ±7 


4.7 ±0.6 ±0.2 






K+tv- 


37.5 ±0.8 


589 ± 30 ± 17 


17.9 ±0.9 ±0.7 


-0.102 ±0.050 ±0.016 


[-0.188, -0.016] 


K + K~ 


36.2 ±0.8 


1±8(< 16) 


< 0.6 
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